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Abstract. We give the spectral representation for a class of selfadjoint discrete graph Laplacians A, with 
A depending on a chosen graph G and a conductance function c defined on the edges of G. We show that 
the spectral representations for A fall in two model classes, (1) tree-graphs with Af-adic branching laws, 
and (2) lattice graphs. We show that the spectral theory of the first class may be computed with the use of 
rank-one perturbations of the real part of the unilateral shift, while the second is analogously built up with 
the use of the bilateral shift. We further analyze the effect on spectra of the conductance function c: How 
the spectral representation of A depends on c. 

Using Ag, we introduce a resistance metric, and we show that it embeds isometrically into an energy 
Hilbert space. We introduce an associated random walk and we calculate return probabilities, and a path 
counting number. 
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1. Introduction 

One fascination with operator theory is its connections to other areas such as geometry and discrete 
analysis: Because of applications to electrical networks, to statistical mechanics, and to fractals (see e.g., 
[Str06l |JP08bl [JP891 IJSWM IPow76l IPow79j ) . there is a recent increased interest in detailed spectral repre- 
sentation for operators on infinite graphs. In addition to the applications, these connections further suggests 
a search for more a direct link between, on the one hand, metric geometry of infinite graphs, and on the 
other, a spectral analysis of associated families of operators. 

However as we see, classical methods break down in infinite discrete models: With Fourier analysis in 
classical potential theory, it is often possible to represent Laplace operators by multiplication, and hence 
realize the spectral representation this way. The analysis then breaks up into the study of discrete and 
continuous parts. However we show that analogues of this that adapt to the discrete case have strong 
limitations. New tools from operator theory are needed: For example in the discrete case, typically there is 
not a natural Fourier duality available, and the graph may not even be associated with a group in a way that 
facilitates computation of spectral representations. Further, in the study of Laplacians on infinite graphs G, 
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there are several Hilbert spaces in the picture. Choices must be made: There are Hilbert space completions 
of functions on the vertices in G, and similarly for functions on the edges in G (Definitions [2l [4~Tjl . the energy 
Hilbert space. 

Here we focus on classes of graphs G that require new tools. Our conclusions (Theorem 13. 26[ Corollary 
I4.6( and Theorem I4.11j) imply that not only is there a direct connection between the spectrum of the graph 
Laplacian and the metric geometry of G; but this connection carries over to the detailed fine-structure 
of the multiplicity configurations for Aq. 

Our proofs rely on a mixture of operator theory (section 3) and complex analysis (section 2). Since we 
address three different audiences, for the convenience of readers, we have included a few details which may 
be known to operator theorists but not to graph theorists, and vice versa. 

In this paper we study the operator theory of infinite graphs G, and especially a natural family of Laplace 
operators directly associated with the graph in question. These operators depend not only on G, but also 
on a chosen positive real valued function c defined on the edges in G. In electrical network models, the 
function c will determine a conductance number for each edge e; the conductance being the reciprocal of the 
resistance between the endpoint vertices in the edge. Specifically if e = (xy) connects vertices x and y in G, 
the number c(e) is the reciprocal of the resistance between x and y. Hence prescribing a conductance leads to 
classes of admissible flows in G determined from Ohm's law and Kirchhoff 's laws of electrical networks. This 
leads to a measure of energy directly associated with the graph Laplacian. There are Hilbert spaces H{G) 
which offer a useful spectral theory, and our main results concern the spectral theory of these operators. In 
a recent paper |Jor08] it was proved that the graph Laplacians are automatically essentially selfadjoint, i.e., 
that the associated operator closures are selfadjoint operators in H{G). 

Here we give a spectral analysis of the graph Laplacians. We are motivated by a pioneering paper |Pow76] 
which in an exciting way applies graphs and resistor networks to a problem in quantum statistical mechanics. 

There are many benefits from having a detailed spectral picture of graph Laplacians: We get a spectral 
representation realization of the graph Laplacians, the operators Aq, c , i-e., a unitarily equivalent form of 
these operators which may arise in a variety of applications. See e.g., }Arv02| [PS72] . 

In Corollarv l4.8( we obtain a candidate for Brownian motion (independent increments) on a general graph, 
and for the tree graphs we show in Proposition ^. 51 that the increments go like the square root of the distance. 

In the course of the proofs of our main results, we are making use of tools from the theory of unbounded 
operators in Hilbert space, especially [Voi85] . but also von Neumanns deficiency indices, operator closure, 
operator domains, operator adjoints; and extensions of Hermitian operators with a dense domain in a fixed 
complex Hilbert space. References for this material include: J0r77l p0r78[ IJP00[ INel69. vN3l| [Sto51j . For 
analysi s on infinite g raphs and on fractals, see e.g., |BHS05l ICS07l lDS8l IDJ061 IHKK021 IHut81l IJP981 
IJKS07l|KTg03llStr06l . 

1.1. Technical details. Let G = (G^ ,G^) be an infinite graph, G^ for vertices, and G^ for edges. 
Every x in G^ is connected to a set nbh(a;) of other vertices by a finite number of edges, but points in 
nbh(x) are different from x\ i.e, we assume that x itself is excluded from nbh(x); i.e., no x in 

G (0) 

can be 

connected to itself with a single edge. Let c be a conductance function defined on G^\ 

Initially, the graph G will not be directed, but when a conductance is fixed, and we study induced current 
flows, then these current flows will give a direction to the edges in G. But the edges in G itself do not come 
with an intrinsic direction. 

A theorem from |Jor08j states that the Laplace operator A = A c is automatically essential selfadjoint. By 
this we mean that A is defined on the dense subspace D (of all the real valued functions on G^ with finite 
support) in the Hilbert space H = H(G) := l 2 (G^). The conclusion is that the closure of the operator A 
is selfadjoint in H , and so in particular that it has a unique spectral resolution, determined by a projection 
valued measure on the Borel subsets the infinite half- line R+; i.e., the spectral measure takes values in the 
projections in the Hilbert space l 2 (G^). We work out the measure. 

In contrast, we note that the corresponding Laplace operator in the continuous case is not essential 
selfadjoint. This can be illustrated for example with A = on the domain D of consisting of all C 

functions on the infinite half-line R + which vanish with their derivatives at the end points. In this case, the 
Hilbert space is L 2 (R + ). 

Our main theorems are as follows: In section [3] we consider a natural graph Laplacian Aq on the infinite 
graph G whose vertices form an A^-fold branching tree. In Theorem I3.16[ we show that Aq has a natural 
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representation in Voiculcscu's Fock space T{Hm) where Hm is an V-dimensional complex Hilbert space. 
The link between Laplacians on trees and the Fock space was first made in |GG05j . 

In Theorems 13.191 13.211 and 13.261 we give an explicit spectral representation of A& with the spectral 
measure being Wigner's semicircle law (see |Wig55| ) or a related measure. Theorem 13.261 accounts for the 
spectral multiplicity of Aq as a selfadjoint operator. In section |4] we introduce the resistance metric for 
graphs in general; and (Proposition 03]) we compute this metric explicitly for Aq corresponding to the case 
when G is the graph of the AT- fold branching tree. In Proposition 13.71 this is contrasted with the analogous 
but much simpler case for the infinite lattice graphs. 

In section^ we establish the connection between the graph Laplacian on one side, and on the other, 
a certain metric d on the vertices G^ of the graph, known in electrical models as the resistance metric. 
Using Ag, we show that the (G^°\d) embeds isometrically into a Hilbert space built directly from Aq, the 
energy Hilbert space. In Proposition 14.51 we calculate this for the case of the tree-graphs. Proposition 14.91 
makes the connection to random walk, return probability, and (Theorem 14. lip to path counting; computed 
from the moments of the spectral measure. 

In a general context graph Laplacians includes discrete Schrodinger operators, see e.g., |ALM07j . and 
Sim95]. The present paper restricts its focus to the graphs G having as vertex set trees built from the 
infinite iteration of an V-fold branching rule; but we contrast our results with those of other infinite graphs, 
for example lattice graphs. Moreover the study of V-fold branching rules is of independent interest, for 
example in signal processing, see e.g., |Jor06 . Another reason for the restriction in focus is that we then are 
able, with the use of Voiculescu's Fock space construction, to write down the complete spectral representation 
for the corresponding graph Laplacian, including a geometric model for the cyclic subspaces occurring in the 
spectral multiplicity table. 

While there is a number of related studies of graph Laplacians in the literature (e.g., |JP08a] L the detailed 
spectral picture has so far received relatively little attention. 

The focus of the present paper is instances of explicit spectral representations. While there is already 
a large literature on graph Laplacians, so far we have only encountered relatively few instances where the 
complete details are worked out for associated spectral representations. Since the geometric possibilities of 
graphs is vast, then so is the associated spectral configurations. Spectral approaches are manifold, and here 
our focus is: formulas for explicit spectral representations. 

1.2. Applications and connections to related results in the literature. A list of recent and past 
papers of relevance includes [Str08| (received after we completed the present first draft); |Car72( [Car 73a) 
ICar73bl ICR061 IChuf)71 ICdV99l ICdVOil lJor83] , and Wigner's original paper on the semicircle law |Wig55| . 
While graph Laplacians have numerous applications (see our cited references) in the theory of fractals 
Str08 , in combinatorics |Chu07j . in random walk models [SC Z07 . in free probability Voi85j, in operator 
theory |HH73j . and in harmonic analysis [RPB071 [Car73bl ICdV99] : they are of relevance as well in math- 
ematical physics; see e.g., [ASWQ61 IAF00I IBre071 IGG051 IMoh911 IRR951 ISoTui] . The mathematical physics 
literature contains a wealth of spectral theoretical results, each with a particular focus. Here we give formulas 
for the spectral transform of the tree-graph Laplacians, including multiplicities. And, conversely we show 
that the tree-graphs reflect themselves directly and in an explicit manner in the spectral multiplicities for 
the associated operator. 

In Voiculescu's original approach to spectral representation, he relied on a high powered result of Helton- 
Howe (traces, index and homology, HH73 ). While one does arrive at a representation this way this argument, 
and the big machinery, leaves one wonder if there is instead a direct and computational way of getting it. 
We find this! 

While authors of earlier more "practical" uses of discrete Laplacians and discrete Schroedinger operators 
have found a variety of spectral features of related operators (see the cited references above) , none of them 
offer a complete spectral transform, and they do not have a spectral duality with Voiculescu's Toeplitz C*- 
algebra, as we do here. We feel that this is of independent interest. It not only extends what was in the 
literature already, but it also offers a unifying framework. 

The latter, plus our explicit formulas, is one of our key points in our analysis and in our derivation 
of complete spectral transforms. In fact we further give a formula for the spectral transform. And we 
use free probability methods from Voiculescu's Toeplitz C*-algebra [Voi85j in pinning down the complete 
spectral multiplicity picture for graph Laplacians on trees. Conversely, we use this in computing moments 
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and resistance distances for the graphs themselves. Thus, our main result shows that the tree graphs reflect 
themselves directly in the spectral multiplicities for the associated graph Laplacian. This may well work in 
both directions: Our spectral theoretic results may be of use in free probability computations. 

2. Definitions 

Definition 2.1. Graph Laplacians. Let G — {G^ Q \G^) be a non-oriented graph with vertices G^ and 
edges G {1 \ If e G G^ we assume that the source s(e) and the terminal vertex t(e) are different. If 
(xy) G G^ we write x ~ y. We further assume that for every x G G^ the set of neighbors 

(2.1) nbh(x) := {y G G (0) | y ~ x} 
is finite. 

Let c : G^ 1 - 1 — * R+ be a fixed function (called conductance). Let l 2 (G^) be the Hilbert space of square 

summable sequences, i.e., 

(2.2) \\v\\f 2 := Y, Wx)\ 2 <oc, 
and 

(2.3) (vi,v 2 ) p := J2 ioYv 1 ,v 2 el 2 (G^). 

Let V be the set of all finitely supported elements in l 2 (G^), i.e., v G V iff there exists F c G' -' ( a 
finite subset) such that v(x) = for all x G G^ -* \ i* 1 . 

Given (G,c), the corresponding graph Laplacian A = Ag jC is 

(2.4) (Av)(s)=$^c(aw)(w(aO-i;(y)), (« G 23) 

It was proved in |Jor08j that A is essentially selfadjoint, i.e., that the operator closure A is selfadjoint. 
Hence for every (G, c), there is a projection valued measure 

(2.5) E A :£([0,oo))->Proj(Z 2 (G(°>)) 

where £> stands for the Borel sigma-algebra, and Proj is the lattice of all selfadjoint projections, i.e., P = 
P* = P 2 . The function E A (-) in (|2.5|) is countably additive on B, and satisfies 

(i) 

E A (A 1 nA 2 ) = E A {A 1 )E A (A 2 ), (A 1 ,A 2 eB); 

(ii) 

Jp = / ^A(rfA); 
Jo 

(iii) 

/» CO 

A= / XE A {d\). 
Jo 

Moreover if / is a Borel function, the operator /(A) is given by functional calculus, 



/(A) = J f(X)E A (dX), 
and a vector v G / 2 (G^ ^) is in the domain of /(A) iff 

(2.6) / |/(A)| 2 ||.EA(<iA)t>|| 2 < °° 

Jo 

For«o G Z 2 (G<°)), ||vo|| - 1, set 

(2.7) MoO) := II-BaOoII 2 = («o , £a(>o> 

. We call fig the spectral measure associated to A and i/ie vector vq. See also Lemma 12.31 below. 
For a measure /iq on R, define 

(2.8) F M0 (z) := [ ^—dfioix), 
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the Borel transform, z G C \ R. 

Lemma 2.2. If Ra(z) '■= (A — zlp)^ 1 is the resolvent operator, then 

F^z) = {v , Ra(z)vo), zeC\R. 

Proof. Using (iii), we have 

Ra(z) = [ {x- z)- 1 Ea(<Ix), 

and for z G C \ R, 

{vq , Ra(z)v )i2 = / (x - z)" 1 (v , E A (dx)v ) = / (x - zy 1 \\EA{dx)v \\ 2 = (x - z)' 1 dfj, (x) = i^ (z). 
JR Jr Jr 

□ 

Lemma 2.3. Lef «o and /io fee as above, and let 7{a(vo) C i 2 (G^ ^) fee i/ie A-cyclic subspace generated by 
Vq. Set 

H A (vo) 3 /(AH K /(•) G L 2 ( Mo ). 

T/ien W extends to a unitary isomorphism of H.a(vq) into L 2 (/zo) which satisfies the intertwining relation 
M X W = WA on Ha(vq). 

Proof. See e.g.. |Nel69ll5to90] . □ 

Definition 2.4. Intertwining operators. 

Let Hi, i = 1,2 be Hilbert spaces, and let W : 7ii — > be a unitary isomorphism of Hi into 7^2- Let 
JFi, z = 1,2 be families of operators in the respective Hilbert spaces. We say that W is intertwining, or 
equivalently that the two families are unitarily equivalent if there is a bijection ip : T\ — > F2 such that 

(2.9) WT = p{T)W, (Te.Fi). 

Definition 2.5. Cuntz relations. 

Let H be a Hilbert space of finite dimension N. Let Ti be a second Hilbert space, infinite dimensional, 
and let B(TL) be the C*-algebra of all bounded linear operators in Ti. The norm in B(TL) will be the operator 
norm, i.e., if T G B(H) 

(2.10) ||T|| :=sup{\\Tv\\\veH,\\v\\ = l}. 
A linear function p : H — ► B(Tt) with the property that 

(2.11) <p(v)*<p(v) = \\v\\ 2 I n , (veH). 

is called a Cuntz- Toeplitz representation, and the operators {<£>(«) | w € i?} are said to satisfy the Cuntz- 
Toeplitz relations. 

Lemma 2.6. Let p : H — ► B(Tt) be as in Definition \2.5\ Then p is a Cuntz- Toeplitz respresentation if and 
only if for all orthonormal basis \i) in H , i — 1, . . . , N , the operators 

(2.12) Zi :=¥>(!*» 
satisfy 

(2.13) T*Tj = and 

Proof. Let v & H, and set w := 5Z i=1 fi | i) - Then by linearity = X)i=i w »Zij see f)2.12[) . Moreover, using 
(|2.11jl and its polarization, we get 

JV AT AT 

(2.14) £ h| 2 = ||^)|| 2 = \\p(v)*<p(v)\\ - EE^* 2 ? T *' ( u G 

i=l 3=1 fe=l 

The relations in (|2. 13[) are directly equivalent to this. 

□ 
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Definition 2.7. Let TV 6 N, TV > 2, and let {Ti}f =l be a family of operators in a Hilbert space 7i. We say 
that this is a representation of the Cuntz algebra On if (|2. 13[) holds and, 



N 



(2.15) J2 T ^ T i =I 

i=i 

We say that this is a Toeplitz system, or a representation of the Toeplitz algebra Tn if (|2.13[) holds, but 

N 

(2.16) Y. T i T i <I 



Definition 2.8. Unilateral shift. 

Let H be a complex Hilbert space and let T : TL — > TL be an isometric operator, i.e., an isometry. We say 
that T is a (unilateral) shift if one (and hence all) of the following equivalent conditions are satisfied: 

(i) There exists a Hilbert space H and a unitary isomorphism W : H — > such that W inter- 
twines T with 

(2.17) (x ,xi,x 2 , ■ • ■ ) -> (0,a;o,xi,ar2,...) 

on £o JT. Then dim if is called the multiplicity of T. 

(ii) For any a; G H, lim^oo T* n x = 0. 

(iii) The projections P n := T n T* n satisfy P x > P 2 > ■ ■ ■ > P n > P„+i > ■ ■■ and inf„ P n = 0. 
See |SNF70| . 

Definition 2.9. Bilateral shift. 

Let 7i be a complex Hilbert space and let T : H — > H be an isometry. We say that T is a bilateral shift 
if one (and hence all) of the following conditions is satisfied: 

(i) There exists a Hilbert space H and a unitary isomorphism W : H — > X^oo* ^ such that W 
intertwines T with 

(Xn)n£l — ► (x„_l)„ e z 

Eoo £D rr 
-°° 

(ii) There is a unitary isomorphism VT : TL — > L 2 (T,H) such that intertwines T with the multipli- 
cation operator M : / — > z/(z) on L 2 (T, if). Here L 2 (T,H) denotes the Hilbert space of vector 
valued functions T — > H, measurable and satisfying J" T ||/(z)|| 2 dv(z) < oo, where v is the Haar 
measure on circle T := {z e C | \z\ = 1}. 



3. Hilbert spaces and a graph model 

We begin with a particular graph, called the tree with TV-fold branching. It is formed out of an alphabet, 
say A of size TV as follows. We first form words Wk = Wk(A) of length k, for each k, with letters chosen from 
A. It is often convenient to take A to be the cyclic group Zjy of order N. The tree V = Vn will be the union 
of the sets Wk for k = 0, 1, ... , such that the k = case corresponds to the empty word. We then form a 
graph Gn with vertices V = Vn consisting of the points in the tree with TV-fold branching. Specifically V 
is the union of the words Wk, with the understanding the Wo is the empty word. The edges in Gn will be 
made up of nearest neighbors in V as follows: If x is in Wo, the set of edges emanating at x consists simply 
of the points in A. If x is in Wk , k > 1, then the edges in Gat consist of the lines connecting two vertices, x 
one of them, and the others resulting from x = (x\X2 ■ ■ ■ Xk) by truncation at the tail end and by addition 
of a letter from A also in the tail end. So the nearest neighbors for the vertex x consist of the TV + 1 vertices 
{x\X2 ■ ■ ■ Xk-i) and (xa) where a is chosen from A; so TV + 1 nearest neighbors in all. 

Starting with this particular tree graph Gn there are two naturally associated Hilbert spaces, the first is 
simply the / 2 -space over Vn, and the second is a Fock space T(H) = T{Hn), Hn ■= 1 2 {2,n) = C N which 
carries additional structure. This additional features of the Hilbert space T{H) will be needed later. 

In Proposition ^. 14l we show that the two Hilbert spaces are naturally isomorphic. While this is known, we 
have included a proof sketch in order to fix our notation. The Fock space J~(H) is used in particle physics in 
the conventional description of (infinite) quantum systems with particles governed by Boltzmann statistics. 
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Definition 3.1. The infinite N-aiy tree is the graph with vertices 

V := {0} U {wi . . . w„ | n > 1, ui, . . . , u n G {1, . . . , N}} , 

and edges given by the relations ~ i, and • • • w n ~ u>i . . . u n i for all n > 1, i, u>i . . . u> n G {1, . . . , AT}. 

Thus V is the set of finite words over the alphabet {1, . . . , N}, including the empty word 0, and we have 
edges between a word u> and the word uii obtained from u) by adjoining a letter at the end. 




Definition 3.2. Define the maps <r, r<, i € {1, . . . , iV} on V: 

(3.1) <t(0) =0, £r(«i...w n )=wi...w B ^i 1 (n>l,ui,...u n e{l,...,N}) 
Thus ct(z) = for alH € {1, . . . , TV}, so a removes the last letter of the word. 

(3.2) Ti($)=i, Ti(w 1 ...w n )=wi...u n i, (wi,...,w n G {l,...,N},i£ {1,...,N}) 

So n adjoins the letter i at the end of the word. 

Define the operators U and Si, (i G {1, . . . , N}) on l 2 (V) by 

(3.3) Uv — v o cr, S{V — v o Ti, (i G {1, . . . , iV}) 

In the following we show that the operators S{ the i 2 (F) as an operator system generate a representation 
of the Toeplitz algebra in several variables (see |Voi85j .) In |Voi85j this Toeplitz algebra is realized in the 
Fock space T(H). The Hilbert space F{H) carries a system of creation operators T v indexed by vectors v 
in H. 

As shown in formula (|3.7ll . the corresponding adjoint operators T* are annihilation operators. In analyzing 
these operators, Dirac's bra-ket notation is convenient. We will select an orthonormal basis (ONB) {\i)} in 
H, and denote the corresponding creation operators Ti , i.e, T{ = Tuy 

In Proposition [3HH we show that the unitary isomorphism W : l 2 (V) — > J-(H) intertwines the operators 
Si with the corresponding system T* in T(H). 

This helps us identify a universal Hilbert space representation for iterated function systems (IFSs), and to 
make precise the quotient of the Toeplitz algebra by the compact operators as a Cuntz algebra (see (Jor04 ) . 

The next two lemmas require some simple computations. 

Lemma 3.3. For f G l 2 (V), i G {1, . . . , N} one has: 
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(i) For alluJi,...,Lu n g {1,...,JV}; 



or n q c J ^lji ...tOn—i 1 ^/ ^ r r 

U, 2/ CJ n 7= X 



(ii) For all CJi, . . . , cj n G {1, . . . , N}: 



(iii) For a// wi, . . . , w n G {1, . . . , N}: 

Si 3<6 = ^ui...u„ — 0~Ui...u>„i- 

Lemma 3.4. The operators Si satisfy the following relations: 

(i) S t S* =1 for alii e {1,...,JV}. 

(ii) Let P% be the projection in I (V) onto the vector 5$. 

N 

i=l 
l2l 



(i) u* = Eti St + p 



Proposition 3.5. Let P$ be the projection in I (V) onto i/ie canonical vector 5$ 
(i) 

A = (iV + 1)1 - (]T 5, + (7) = (N + 1)7 - (1/ + 17* - P t ) = (N + 1)7 - (£ S. t + £ 5* + P ). 



JV TV JV 



= 1 i=l i=l 

2( 



Proof, (i) For v, v' in £ (V), we have 

{Uv,vf)=v(H)i/(9)+ v(u 1 ...u n -i)v'(uii...uj n ) = 



n>l,Wi ,. . . ,LOn 

«(0V(0)+r" ^ ' 



2^u'(wi)+ ^ i>(wi...w n _i)^u(wi...w n -iw tl ). 

wi oji . . . ..u;™ — i 

Therefore 

N N 

(U*v)(<b) = «(0) + ^ «(i), . . . O = XI «(wi ■ • ■ 

it C/* 

(ii) We have 



»=i i=i 



This implies that U* = J^fLi S i + Pt 



N N 



(A«)(0) = JVu(0) - ^ «(i) = (JV + 1)«(0) - (^ «(*) + «( 
and 

JV 

(Au)(wi . . . w n ) = (TV + l)v(wi . . . w„) - (v(lji . . . £J n -i) + X w ( Wl ■ ' ■ 



This implies that A = (JV + 1)7 - Si + U). Tne other equalities follow from (i). □ 

Remark 3.6. In equation (|2.4p we introduced our general class of graph Laplacians Ag, c . As noted there, 
there is a Laplace operator for each graph G and for each choice of conductance function c. Our paper is 
about spectral theory of the graph Laplacians, and the spectrum depends on both the graph, and the choice 
of conductance function. 

To understand the graph dependence, we may contrast two cases: Case 1: the JV-fold tree graphs (Defi- 
nition [33] and Figure [1]) and the use of the unilateral shift (Definition 12. 8|l . Case 2: Lattice graphs (details 
below). For the latter we show that the spectrum is determined by the bilateral shift (Definition 12.91 ) 

Let G = (G (0) , G (1) ) be the rank d group with vertices G*- - 1 = Z d , and edges defined by n = (ni, . . . , n<j) ~ 
m = (mi, . . . , rid) iff there exists k € {1, 2, . . . , d} such that \m,k — n&| = 1 and mj = nj for j =/= k. 
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For v G l 2 (Z d ) set 

(3.4) (Av)(n) := ^ v{n) ~ v(m); 

for example, if d = 1, 

(At>)(n) = 2v(n) - u(n - 1) - + 1). 
These operators are used in numerical analysis, in electrical network analysis, in electrical network models 
in physics; (see e.g., [Kig03, Pow76|1 but with a variety of choices of the conductance function; see (|2.4[) 

Proposition 3.7. Let d G N. The graph Laplacian A in (|3.4p has the form I — 2ReT ( j where KeTd = 
\{Td + TJ), Trf = <g)fT and T is the bilateral shift; see Definition \2.9\ 

Proof. Introducing the Fourier transform on T d (the d-torus), we get the unitary equivalence 

(3.5) V(z) = V nZ n «-> (Vn)n G l 2 ^) 

nez d 

and by Parseval's equality 

IMIz,2(TN) = l U "| 2 - 

Substituting (|3.5|) into (|3.4p . we get 

d d d 

(Av)(z) = (27V-^z fe -^z fc )«(z) = (2JV -£)2Re 

fc=l k=l k=l 

Or setting 0fe = e lXk , fe = 1, . . . , d, 

d d 

(Av)(x) =2Y J {^-cosx k )v{x) =4^sin 2 (^\ v (x). 
fc=i fe=i 

□ 

Remark 3.8. There are several contrasts between this case and the Laplace operators of tree graphs, see 
Theorem 13.261 One is that the spectral measure here simply is a pull-back of the Haar measure on T d with 
the function given above; second that there is no rank one perturbation; and third that the spectrum is 
simple in the present bilateral case. 

We also note that the Laplace operator for the lattice case has simple spectrum while the case of the tree 
has an intricate multiplicity structure, see Theorem 13 . 1 61 for specifics. 

3.1. The Fock space. 

Definition 3.9. Let H := C N and be a fixed unit vector. Then the Fock space is the Hilbert space 
T:=T{H) := Cil @ H @ (H ® H) ® ■ ■ ■ @ H g ■ ■ ■ ® H . . . . 

n times 

For a vector v G H we define the operator T v on T by 

(3.6) T v x := x ® v, x G T, {T v il := v). 
Proposition 3.10. For v G H, the adjoint of the operator T v is given by the formula 

(3.7) T„*rj = 0, T*{xi ® ■■■®x n ) =xi O ••■®x n - 1 {v, x n ) , (n > 1, x x , . . . ,x n G H), T*n = 0. 

The following Proposition summarizes a number of geometric properties of our isometries in the Fock 
space. We begin with the operators T and T* for T — T v . In Definition 12.71 we emphasize that a unilateral 
shift is really an isomorphism class; referring to unitary equivalence, and that the multiplicity is a complete 
isomorphism invariant. Specifically, in the Proposition below, we spell out a particular shift representation 
for the operator T v when ||u|| = 1; and this refers to the Fock space which we will need later. In fact, for a 
fixed v, we must identify the closed subspace in the Fock space which is shifted by powers of the isometry 

We use Voiculescu's framework |Voi85| and [OraOl . Moreover we include here the details we will be using 
later in the proof of our main conclusions regarding spectral representations of graph Laplacians;- in the case 
of trees, these are operators derived from the T v system. 
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Proposition 3.11. Let N G N, and let H be a Hilbert space of dimension N. Let T := F{H) be the Fock 
space (see Definition \3.9[ ) Let v G H satisfy \\v\\ = 1. Then T = T v is a shift operator in T . The multiplicity 
of T is as follows: mult(T) = 1 iff N = 1. And if N > 1, then mult(T) = oo. 

Proof. It follows from Proposition 13. 101 that T is isometric, T : T — ► and that 

ne N{T*) = {£eT\T*£ = Q}. 

If N = 1, then N(T*) = Cf2; and otherwise N(T*) contains in addition the following vectors in T: £i ® 
•••(g) Cfe-i ® ik where & G ff, and £k -L «, i-e., (6, , w ) = 0. 

Moreover, together with fi, the closed span of these vectors exhausts N(T*). 

It is true in general for a fixed isometry (T, Tt) that the direct sum Hilbert space 



(3.8) J2 N(T*) 



k=0 

is isomorphic to a closed subspace in Ti. Specifically, if x k £ N(T*), then 

oo 

(3.9) ^T k x k -^xq®x 1 @x 2 ® ... 

fc=0 

is a well defined isometry, i.e., it satisfies 

oo oo oo © 

(3.10) iiE Tfc ^n« = EiNii 2 = iiE ^n 2 <^ 

fc=0 fe=0 fe=0 

This is a consequence of the following identities: 

(3.11) (TS Xj , T k x k ) H = 5 jtk ( Xj , x k ) , (j, k G No, ay G N(T*)). 

Introducing P k := T*T* fc from (pT9> ([330)1 . we see that a vector £ G H satisfies inf fe \\P k £\\ 2 = iff f G 
V fe Ran(T fe ) = (A fe jV(T* fc )) x . 

Set H = T = T{H), and assume AT > 1. Set T — T v , for v £ H, \\v\\ = 1. Then 

(3.12) Pfc+itf^ =0, (fee N ) 

with the convention i?® = CO. 
If k < m, then 

(3.13) PkU^H = {I®™-*) ®\ v®- - -®v) ( v®- - -®v \ 

k times k times 
where we use Dirac's notation |-)(-| is denoting rank one operators. 
Since T = J^kxj^ H Xk > we nave the following representation: 

fc>0 

and 

(3-14) ll£lll=^IIU|™H 

m>0 

Using ([3"32"j) - ([3TT3) . we then get 

ll^eil 2 = ]T W P ^n\\ 2 < H^ll 2 ^ as fc -> co; 

m>k rn>k 

where we used (|3.10[) in the last step. 

Hence T is a unilateral shift according to Definition 12. Sf iii). 
It follows from (|3"T3"|) that mult(T) = 1 iff N = dimH = 1. 

If N > 1, then N(T*) is the closure of the span of and the infinite sequence of closed subspaces 

(®qH) ®(He {«}). 
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Hence this sum is an infinite dimensional Hilbcrt space. We used the terminology: 

H 6 {v} := {w e H I (w , v) = 0}. 

□ 

Definition 3.12. Let H be an TV-dimensional Hilbert space and let T = J~(H) be the Fock space. Let 
v G H, \\v\\ = 1 be given. Set Ran(T l) ) =the range of the isometry T v . Then T V T* is the projection onto 
Ran(T„), and / - T V T* is the projection onto N(T*). 

Corollary 3.13. Let H and T = T{H) be as in Provosition \3. 1 1[ Suppose N — dhnH > 1, and let v € H, 
\\v\\ = 1 be given. Then 

(3.15) Cft® Ran(T w ) = N(T*) 

w£HQ{v} 

Proof. Set v\ — u, and extend to an ONB: V2, ■ ■ . , vn for H . Then we have 

JV 

(3-16) I? - T^r:, =Pn + Y, T ^K 

i=2 

But we noted that I? — T Vl T* t is the projection onto N(T Vx ), while the right-hand side in (|3. 16[) is the 
projection onto 

w±.v 

□ 

Proposition 3.14. Consider the infinite N-ary tree in DeHnition \S.l[ Let (ei)fL 1 be an orthonormal basis 



for H — C . And let T = T(iT) be the corresponding Fock space. The map W : l 2 (V) — > T defined by 

(3.17) W(<y ) = n, yV{8 Ul .. Mn ) = e Ul »'--®e Un , (w u ...,u n e{l,...,N}) 

defines an isometric isomorphism between the two Hilbert spaces. 
Let sa '■= ^2iLi e i- Then 

(3.18) WSiW*=T*., (ie{l,...,JV}), 
and 

(3.19) WAW* = (N + l)L - (T SA + T S * A + P u ) 



where Pq is the projection in T onto the vector f2. 

Proof. The map W defines an isometric isomorphism because it maps the orthonormal basis into 
an orthonormal basis {SI} U {e ull ® • • • ® e Un \ wi,...,w n £ {1, . . . , N}, n e N}. 
Next, we check equation (|3.18|) on this orthonormal basis. 

WSiW*Q = WSi5 $ = W5 $ o n = o = r*o. 
WSiW*e Wl O • • • <g) e Wn = WSj^...^,, = W^...^ o t, = W Wl ... tl) „_ 1 (5 w „, i = 
e Ul <8 • • • ® e w „_ 1 (e Wn , e 4 ) = r*.e Wl (g) • • • (g> e Wn . 

This implies (j3~T5| 

It is clear that WP®W* = P n . 

From equation (|3. 1 8[) and Proposition 13 . 51 it follows that 

N N 

WAW* = (N + 1)1 - ($2 T*. +^T ez + P u ). 

i=l i=l 

But since T ei = T SA , equation (|3~Tg|) follows. □ 
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3.2. Decomposition into cyclic subspaces. 

Definition 3.15. Let {xi, . . . , xjv} be an orthonormal basis for H = C N with the first vector X\ = Sq := 
I I^H = -^sa- (sa is defined in Proposition 13. 14|) . Define the following subspaces of T: 

Hn •= span{£7, (g p k=1 s Q \p> 1}, 

Hi!...,, := spmix^ (g> ■ ■ ■ <g> x iri (g> ® p k=1 s \ p > 0}, ii,. ..,i n G {1,. . . ,N},i n ^ l,n > 1. 
so Xj n should not be the vector in the basis corresponding to sa- 

Theorem 3.16. XTie cyclic subspace decomposition of the graph Laplacian. 

(i) TTie subspaces Tin and Ti,i 1 ...i n , ix, ■ ■ . , i n G {1, ■ ■ • , N}, i n =/= 1 are mutually orthogonal, 

Tin © ®ii,...,i n e{l,...,N},i n ^l7~ii 1 ...i n = J 7 , 

and they are cyclic subspaces for the operator WAW* . 

(ii) Let S be the unilateral shift of multiplicity one, 

S : l 2 (N ) -> l 2 (N ),S(x ,x 1 ,x 2 ,.. .) = (0, x , x t ,x 2 , ■ ■ ■ )• 
TTie restriction of WAW* to 7Yji is unitarily eguivalent to the operator 

(3.20) £> n := (JV + l)/p (No ) - 2\/iVRe 5 - Pa . 

For all i\, . . . , i n S {1, . . . , AT}, i„ 7^ 1, and n > 1, £/ie restriction of WAW* to Ti^...^ is unitarily 
equivalent to the operator 

(3.21) Z? := (JV + 1)J; 2 (f}) — 2\/]VRe 5. 

Where KeS = S+ 2 S , and P$ is the projection in I 2 (No) onto the vector Sq. 

Proof, (i) To prove that the spaces are mutually orthogonal, take x^ ® • • • ® Xi re ® ®^ =1 So an d sCjj ® • • • <8> 
Xj m (g) g> q k=1 so. If the lengths are different, i.e. n + p 7^ m + q, then the two vectors are orthogonal. If 
n + p = m + q and n < to, then since j m 7^ l, it follows that so -L Xj m , and again the vectors are orthogonal. 
Similarly, if n > to. If, in addition n — m, then if i\ .. .i n 7^ ji . . .j' m then ik =/= jk and Xi k _L Xj fc and the 
two vectors are orthogonal. 

To prove that these subspaces span the entire space, we have that {f2, x^ ® • • • ® Xi n | i& G {1, . . . , iV}} is 
an orthonormal basis for T . If all ik = l, then Xj x £5 ■ • • (8> Xi„ is in 7in. If some 7^ 1, take the last one as 
such, and 2^ ® • • ■ <8 Xi n G T~U 1 ...i h - 

It remains to prove that these subspaces are cyclic for WAW*, or equivalently by Proposition 13.141 for 
A := T SA + T S * A + P Q . Note that T SA = VNT So , so A = yfN(T Sa + T S *J + P n . 

Take n > 1. Take Xj t ® • • • £3 Xj n with ii, . . . , j„ G {l, . . . , ./V} and « ra 7^ 1. Then 



A(xi 1 ® • • ■ ® sci n ) = v ATxij ® • • • ® Xi n Cg) so + + 0. 

Assume by induction on p that x^ (8 ■ ■ • <8> x^ n ® (S>? =1 So is in the cyclic subspace of A generated by the 
vector Xi x ® • • • ® Xi n . Then 



Aa;^ <g> • • • (g> Xj„ ® <8fe =1 so = v iVxi! ® • • • ® Xi n ® <8>f =1 so + v iVxij <8) • • • ® Xj n (g) <8>fc =1 so + 0, 

so Xjj ® • • • ® Xi re <g> <8)^}so is in the same cyclic subspace. This implies that Hi 1 ...i n is the cyclic subspace 
of the operator A (hence of WAW*), generated by the vector x^ <g> • • • <g> Xi n . 
Similarly for Tin- 

(ii) For i\ . . . i n G {1, . . . , N}, i„ 7^ 1, define the operator Wi 1 ...i 7l : / 2 (No) — > 'H% 1 ...i n i 

yVi 1 ...i n 5 P = x\®---®x in ® ®^ =1 s , (p > 0). 

The operator W^...!,, is clearly unitary since it maps an ONB to an ONB. A simple calculation shows that 
"WiL.An intertwines D and W*AW on Ti.i 1 ...i n - Similarly for Ha- 

□ 
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Lemma 3.17. [Voi85 The spectral measure o/ReS*. Let S be the unilateral shift of multiplicity one on 
/ 2 (No). Let /x c be Wigner's semicircular measure on [— 1, 1], i.e. [i c is absolutely continuous with respect to 
the Lebesgue measure on [—1,1], and 



d/ic = ( fVl - X s , ifxe [-1,1] 
dx \ 0, otherwise. 

There exists an isometric isomorphism $ : / 2 (No) — * £ 2 (/i c ) 7 such that $(<$o) = 1 (the constant function 
1 ), and 

$ReS" = M x $, 

where M x is the operator of multiplication by the identity function x on L 2 (fi c ), (M x f){x) = xf{x). 

Proof. The proof in [Voi85 involves some heavy machinery (the Helton-Howe formula [HH73]), so we will 
give here a more elementary proof. 

First define the isometric isomorphism J\ : ^ 2 (No) — ► L 2 ([0, ir], — dx), 

(Ji5 k )(x) := sin(fc + l)x, (x G [0, tt], k > 0). 

Elementary Fourier theory (for odd functions on [— tt,tt]) shows that this is an isometric isomorphism. We 
claim that 

(3.22) JiReS = M cosx Ji 

where M cosx is the operator of multiplication by cosx on L 2 ([0, tt]). We have 

JiKeSSo = Ji~^$i — 2 sm (2x) — cos x sin x = M cos x J\ 6 Q . 



For k > 0: 



Ji Re S6k — J\-^{Sk+i + <5fc-i) = 2"( sm (^ + tyx + sinfcx) = sin(fc + l)xcosx = M cosx Ji5k- 



Since J\5q — sinx, we define J2 : L 2 ([0, n]) — » L 2 ([0, ir], ^ sin 2 xdx) 



(J2f)(x) = f(x)^—, (xG[0,tt]). 
sin x 



□ 



Then J2(Ji<5o) = 1, the constant function 1. And clearly M cosx J2 — JiM cosx , and Ji is an isometric 
isomorphism. 

Finally, define the change of variable operator J3 : L 2 ([0, n], ^ sin 2 xdx) — > L 2 ([— 1, 1], /u c ) 

(J3/)(^) = /(cos- 1 a ; ), (are [-1,1]). 

Then J^M cosx = M X J^. To check that J3 is an isometry, use the change of variable cos -1 x = y, i.e. 
x = cos y: 

F— 1 _1 /-l 

/(y) sin 2 y = / /(cos -1 x)(l — x 2 ) -. dx = / /(cos -1 2))Vl — a; 2 dx. 

Ji VI — x 2 J-i 

J3 is also bijective, because its inverse can be explicitly computed: change back the variable. 
The desired isometric isomorphism is $ := J3 o J 2 o J-y, 

Definition 3.18. Let /x c the semicircle measure as in Lemma f3. 171 Define the following operators on L 2 (/x c ): 

E(f) = [ f(x)dfi c (x) = — f f(x)yfl — x 2 dx, ( considered as a constant function ) (/ G L 2 (/i c )). 
J-i Ti-J-i 

(A A /)(^) = 2ViVx/(x) + E(f) (x G [-1, 1], / G L 2 (^c). 

(M x f)(x) = xf(x), (xe [-i,i],/ei 2 y. 
Theorem 3.19. i*br N >2, the Laplacian operator A is unitarily equivalent to the operator 

{{N + 1)7 - A A ) © ®« G N((A r + 1)/ - 2\fNM x ), on the Hilbert space ® n > i 2 (^ c ). 
.For AT = 1, £/ie Laplacian operator A is unitarily equivalent to 

(N+1)I-A A on L 2 {ii c ). 



14 DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 

Proof. The conclusion follows directly from Theorem 13.161 and Lemma 13.171 The only thing that remains 
to be proved is that the projection P$ in i 2 (No) is mapped onto the operator E on L 2 (jj, c ) by the unitary 
equivalence in Lemma 13.171 But this is clear since this unitary maps into the constant function 1. □ 



Proposition 3.20. The spectrum of the Laplacian A is [N + 1 — 2\/N,N + 1 + 2yN]. The Laplacian A 
has no eigenvalues. 



Proof. By Theorem 13 . 191 it is enough to analyze the operators and 2\JN M x as in Definition 13.181 The 
spectrum of 2\jNM x is [—2vN,2yN]. Since the measure /i c has no atoms, the operator 2\f~NM x has no 
eigenvalues. 

We turn now to the operator Aa = 2yNM x + E. Since is selfadjoint, the spectrum is contained in 

R. 

We will need the following 

(3.23) l x =- [ l —^du, c ±\, if |A| > 2\fN, Ael. 

7T J -i A - 2vNx 



To prove (|3.23|) . note that the function A i— ► I\ is decreasing on the intervals (— oo, —2>/N] and [2y/~N, oo). 
Also, by a change of variable x — —y we obtain that I-\ = 
Therefore it is enough to compute 

2 f 1 1 — x 2 If 1 ll + x ll + x 

In ntr = — / — 1= dx = == / \ olx = ( use the substitutioni / = u) 

2 ^ nJ-i2VN-2VNx ttVnJ-i\1-x V V 1 - x ' 

4 f°° u 2 4 r 12 . 2 1 

= — = / — = rr du = ( use the substitution u — tan t) = — = / sin tdt = — = < 1 . 

ny/N J (u 2 + l) 2 iry/N Jo ^/N 

This implies (j3~23)) . 

We prove that for |A| > 2\/N the operator XI — A& has a bounded inverse. Let g G L 2 {p c ). We want to 
solve 

Xf(x) - 2VNxf(x) - E(f) = g(x), (x G [-1, 1]). 

Equivalently 

(3.24) f(x)(\-2s/Nx)-E(f)=g(x) 
Integrating with respect to /i c we obtain 

9{x) 



A - 2y/Nx 



E(f)(l-h) = £ 

With dug]) 

(3.25) E(f) = — ^—E I 9 ^ 



dfj. c (x). 



1 - h \ A - 2\JNx 
Then, from (|3~23|) and (|3~2"5)) we get 

(3.26) m = m±m t 

Since |A| > 2y/N the function - 2 ]^ x is bounded on [—1, 1], Therefore, from (I3.25[) . using Holder's inequality 

IW)|| 2 < C\\g\\ 2 , with C depending only on A. Then, from §M, Wfh < C(\\gh + < C'hh, 

with C" depending only on A. This shows that the operator XI — A& has a bounded inverse. 
Next we take A G (— 2yfN, 2y/~N) and we show that the operator XI — Aa is not onto. Take 

9i) ~\ if l--^l^i 

Clearly g is in L 2 (p c ). Suppose there exists / G L 2 (/i c ) such that Xf — A&f — g. Then, as in (|3.26p . 

(3.27) / 9 -^= E ^ f l- , (ze[-l,l]). 

V ' X-2VNx X-2^Nx 
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But the two functions on the left are in L 2 (/i c ) (since g is zero around the singularity x = \/(2yN)), while 
the one on the right is not, unless E(f) = 0. But if E(f) = then (|3 . 2T[) implies that f(x) — around 
x = A/(2v / 7V) and f(x) = 1 otherwise. Then E(f) > 0, a contradiction. 

Thus the spectrum contains {—2yN, 2yN) and since it is closed, it follows that it is equal to [—2y/N, 2yN]. 

Also, if A is an eigenvalue for then as we have seen above, |A| < 2\/~N. If |A| < 2>/~N, and / is an 
eigenvector, then 



But the function on the right is not /^-square integrable (at x = X/{2y/N)), unless E(f) = 0, in which case 
/ = 0. Thus there are no eigenvalues. 

□ 

3.3. Rank-one perturbations. Since our graph considerations are global in nature, on the face of things, 
it may seem surprising that the spectral theory of associated graph Laplacians will involve rank-one per- 
turbations of selfadjoint operators; usually thought of as "local". Similarly, one might think that such 
perturbations might be "harmless", but nonetheless they can be drastic from a spectral theoretic view- 
point; and they are intimately tied in with such deep theories as Krein's spectral shift formula, Aronszajn- 
Donoghue's theory for rank-one perturbations, and the generation of singular continuous spectra. Rank-one 
perturbations further explain resonances and spectral-shift near Landau levels in atomic physics. They 
can generate large point spectrum, as well as singular continuous spectrum, as measured by a computa- 
tion of Hausdorff dimensions. The reader may get a sense of these intricacies from following references 
[AKK041 LAKK051 1 ALM071 IBBR071 IDSS071 IK W021 lPo!98l ldRJLS96l ldRKS06j . 

Theorem 3.21. Define the measure [A c +p on [—1, 1] by 



7T v 

1 - 2N~ 1 / 2 x + N- 

The restriction of the operator W AW* to Ti$i ( see Definition \3.15\ Proposition \3.1J\ and Theorem \3.16\) 
is unitarily equivalent to the operator of multiplication by N + 1 — 2^/Nx on L 2 (fi c+p ). 

Remark 3.22. Our plan is to make use of the moments of the operator Re S from Lemma f3. 171 This refers 
to the cyclic vector introduced in the lemma. As we will see, (|3.30p below, up to a geometric factor the 
moment numbers are the Catalan numbers. Since the generating function for the Catalan numbers is known, 
we are able to use this in analyzing the spectral picture for our particular rank-one perturbations of Re S. 
Details: We begin with a fundamental principle (Theorem 13. 23 j) in spectral theory. The idea behind this is 
based on a fundamental idea dating back to Marshall Stone [Sto90 , and valid generally for the spectrum 
of selfadjoint operators A in Hilbcrt space: One is interested in the spectral measure of a given selfadjoint 
operator A, and its support (= the spectrum of A.) Stone |Sto90j suggested (see Lemma f2 .2[) that spectral 
data may be computed from the study of the resolvent operator R{z) — (A — z) 1 °f A, thinking here of R(z) 
as an operator valued analytic function. Since A is selfadjoint, R(z) is well defined as a bounded operator 
function defined in the union of the upper and the lower halfplane. It is analytic in the two halfplanes, but 
fails to continue analytically across the real axis precisely at the spectrum of A. Evaluation of R(z) in a state 
yields the Borel transform F(z) of the corresponding spectral measure (Lemma 12. 2p . and the same analytic 
continuation/reflection principle applies, hence Theorem 13.231 

As we make a continuation in z across a point x on the real line, Stone suggested (see also Theorem I3.23j) 
that the singular points of the continuation of F(-) occur precisely when x is in the spectrum of A. In the 
absolutely continuous part of the spectrum, the limit from the upper halfplane of the imaginary part of F is 
the Radon Nikodym derivative of the spectral measure. 

This idea was further developed in mathematical physics under the name "edge of the wedge", see e.g., 
[Rud71j . 



Proof. ( of Theorem VJ .21)) By Theorem 13. 161 the operator is unitarily equivalent to (N +1)1 — 2\A/V (Re S + 
2 J r ^ P^ ). We know the spectral picture for the operator Re S. It is given by the semicircular law in Lemma 
13.171 We are dealing here with a rank one perturbation of this operator. We follow the ideas from [Slm95 . 
The moments of the semicircular measure are given by the Catalan numbers (see |RPB07j ). 
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/ 2n 

Let C n := ( _ ) , be the Catalan numbers. 



Then their generating function is (see |Erd06] ) : 



(3.29) C(x) := £ - 1 V 2 [ AX 

n>0 

By the ratio test the series is convergent if < 4. 

The relation between the Catalan numbers and the moments of /z c is (see |RPB07j ): 

(3.30) J ^ x 2n dn c = ^J x 2, Wl~x 2 dx = ^C n 

As follows from e.g., [GLS071 IRPB07j the stated moment relations for the Catalan numbers may be 
derived by induction from the following recursive relations 

k 

(3.31) Co = 1, Ck+i = ^ C n Ck-n- 



In particular 



_ 2(2fc+l) „ 4 fc / 2fc \ (2k 

Ck+1 - k + 2 Cfe ' Cfe -p7v^' Cfe -( fc J-lfe-l 



and 

(3.32) j\™J_^/ R 2_ x 2 dx= Q Ck _ 

Note that the odd moments of the semicircle law are all zero. Formula (|3.29[) follows from (|3.31[) by the 
following simple derivation 

oo oo oo k 

C{x) = C n x n = 1 + xY, C k+1 x k = ( by §M ) = 1 + X J2J2 °nCk~nX k = 1 + xC(x) 2 . 

n=0 fc=0 k=0 n=0 

Solving the quadratic equation for C(x) 1 then yields the two solutions 



1 ± VI - 4x 

We pick the solution from the "-" choice as it satisfies the correct boundary condition C(x) = 1 at x — 0. 



The Borel transform of a measure is 



F(z)= I —dti(x) 
x — z 



We compute the Borel transform of the semicircular measure /i c : 

f 1 —1 f 1 1 —1 1 f 1 

F ( Z )=IZ — : d t J -c( x ) = —l i ^dfx c (x) = — x n dfi c 



x — z z I i 1 — - z ^— ' z 

1 n=0 



— V— — c =— Vc ( — 

n=0 n=0 v 



1 



Thus 



(3-33) F{z) = ^c(±^=-2 Z [l-^-f^ 

For the moment, we know that the relation holds for \z\ > 1, but we will show it actually holds for z € 
C\[-l,l]. 

Theorem 3.23. |Sim95[ Theorem 1.6] Let F(z) be the Borel transform of a measure /i, obeying J ijm < oo. 
Then 
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(i) The singular part of fi is supported on the set 

{x I limlm Fix + ie) = oo}. 

(ii) The absolutely continuous part of \i has Radon- Nikodym derivative: 

^ ac = — limlm Fix + ie). 

dx 7T e|0 

Let F a be the Borel transform of the spectral measure \i a with respect to the cyclic vector 1, of the 
rank-one perturbation M x + aE on L 2 (fi c ). (For us a = See Theorem 13.191 ) 

The Aronszajn-Krein formula is (see |Sim 95, Equation (1.13)]): 

(3.34) F a (z)= F[:) 



1 + aF(z) 



Thus we can explicitly compute F a (z). 
Then 



ImF Q (,)=In 1 F(z)(1 + aF(z)) - lmF{z) 



|l + ai^)| 2 \l + aF(z)\ 2 ' 

By Theorem 13.231 we should have lining — lmF(x + ie) = Wigner's semicircular distribution on [—1,1], 
and outside [—1,1]. 

And it remains only to compute lim e |o |1 + aF(x + ie)\ 2 . Assuming that we get finite limits, this should 
be the distribution of our spectral measure of the perturbation. 

Consider the branch of the square root defined on C \ {x S R | x < 0}. The square root function is 
holomorphic in this domain. 

Take now z € C \ [-1, 1]. Then z 2 e C \ [0, 1], so jg € C \ [1, oo) so 1 - 4y e C \ (-oo, 0]. Therefore the 
function in the righthand side of (|3 . 33[) is analytic in C \ [—1, 1]. Also F is analytic in this domain, therefore 
the formula (|3.33|) is valid for all z € C \ [— 1, 1] . 

The square root: 

As we mentioned above, we define the square root function on C \ (— oo, 0]. An easy computation shows 
that 



1, 


y 


> 


0, 


y 


= 


-1, 


y 


< 0. 



/ — 7— U/x 2 +y 2 + x . y/x 2 +y 2 - x 

(3.35) V x + i y = v- ^ + ts sn(y)v- ^ 

where 

sgn(y) 

' 1 1 1 ( ■ next Lemma can be checked immediately. 

Lemma 3.24. If x < and z n — > x then: 

(i) Iflrazn > for all n then ^J~z~^ — * jy ' —x. 

(ii) Iflm.z„ < for all n then — > —i^j—x. 

Lemma 3.25. If x € [—1, 1], z n = x n + iy n — > x, and y n > for all n, then: 



(i) Ifx> 0, then yjl-jt -> 

(ii) Ifx< 0, then Jl-^r -» -h/^t-1 



Proof. If x > then x n > for n big, and since y n > it follows that Imz 2 > 0. Then Im < so 
Im(l — p-) > 0. Then (i) follows from Lemma f3. 241 
(ii) can be obtained similarly. 

□ 



18 



DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 



With Lemma for x > 0, 



lim F(x + iy) = -2x ( ? ) = -2x + 2i\fi~ 



y|0 



For x < 



limF(x + iy) = -2x ( 1 + i\j ^ ) = -2x + 2i\f\ 



Then, for x G [-1,1], 



.1 , !VT^ |vT 

lim — Im F a (x + iy) 



vio tt v ' (1 - 2ax) 2 + 4a 2 (l - a; 2 ) 1 - Aax + Aa 2 ' 

Then we plug in a = and we obtain that M x + -^=E is unitarily equivalent to multiplication by x 
on L 2 (/j, c+p ) : and the conclusion of the Theorem follows. □ 

Theorem 3.26. Let /j, c be the semicircular measure on [— 1, 1] 

2 

d[i c = — yl — x 2 dx, 

TT 

and let fi c +p be the measure on [—1, 1] given by 



dlLr — - — tit: r dx. 

P 1 - 2N~ 1 / 2 x + N- 1 

Then 

(i) If N = 1, then the Laplacian A is unitarily equivalent to the operator of multiplication by2 — 2x on 
L 2 ([ic+p)- 

(ii) IF N > 2, then the Laplacian A is unitarily equivalent to the multiplication operator 

M c+P 8 ®n=i M c on L 2 {nc +P ) ® ®^ =1 L 2 (^ C ), 

where M c+p is the operator of multiplication by N + 1 — 2^/Nx on L 2 {ii c+p ), and M c is the operator 
of multiplication by N + 1 — 2yNx on L 2 {jj, c ). 

Proof. The Theorem follows directly from Theorem 13. 16[ Lemma T3.171 and Theorem 13.211 □ 



Remark 3.27. Theorem 13.261 shows that the particular rank-one perturbations used in the spectral repre- 
sentation for Aq do not introduce point spectrum. Here G refers to the iV-ary tree. This could not have 
been predicted by the general theory of rank-one perturbations, see e.g., [AKK04] . Nonetheless, as we show 
below, Aq has infinite-energy (i.e., not I 2 -) eigenvectors. For the sake of simplicity, we complete the details 
only in special case of N — 1, but an analogous construction works in general. 

While the Laplace operator Aq for G —the iV-bifurcation graph has absolutely continuous spectrum, the 
following example (for N = 1) shows that Aq may have infinite-energy eigenvectors. 

Example 3.28. On sequences u — (uq, u\, U2, ■ ■ ■ ), set 

(3.36) (Au)o = t«o — tii, an d (Au)„ = 2u n — ti n _i — w n +i, n > 1. 

In general A = is always an eigenvalue. Indeed the vector v = (1, 1, 1, . . . ) satisfies At; = 0. 

The idea in the algorithm of this example is that we generate a finite system of eigenvalues and eigenvectors 
for each n. Given n, the admissible eigenvalues A occur as roots in a polynomial p n . 

Details: Fix n and examine a vector which begins with a finite word w, all letters in w assumed nonzero. 
Letting A be free, and setting v n — 0, you get a polynomial equation v n = p n (A) = 0. Now solve for A, 
and then use the recursion to generate the rest of the coordinates in a A-eigenvector, i.e., use the recursion 
to compute Vk for k = n + l,n + 2, . . . . In each case, when A is fixed, we will get a periodic sequence for 
the entire vector v — (v 0l vi, . . .). Or course this periodicity implies boundedness of (ffc), as k varies in N. 
For each of the admissible values of A you get a one-dimensional eigenspace. So for convenience, we can set 

Vq = 1. 

In the special case of N = 1, set £ = (1, 0, — 1, — 1, 0, 1) and v = (£,,£,,£,,■■■), repetition of the finite word 
£, satisfies Av = v, so A = 1 is also an infinite energy eigenvector. 
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Continuing this process, we can show that for every n€N there is a monic polynomial p n (X) such that 
each root A in p n (X) — is an infinite-energy eigenvalue of A. Moreover the sets of roots are disjoint. 
Specifically, 

pi(A) = l-A, pa(A) = 1-3A + A 2 

and 

Pn+iW = (2 - X)p n (X) -p n -i(X). 

For the case n = 2, the two roots are A± = and we proceed with the analysis of this case. Then 

the eigenvalue problem 

(3.37) Av = Xv 

has non-zero solutions v ^ for the two Golden ration numbers 

(3.38) A = ^ 

In each case, the eigenspace is spanned by the following two periodic sequences There are words £± and rj±, 
£± of length 2 and r\± of length 8 respectively, such that the corresponding A± eigenspaces are spanned by 

(3.39) (£?7, . . . ) infinite repetition. 

Proof. A computation shows that every solution (voViV^ ■ ■ ■) to (|3.37[) satisfies Vi = (1 — X)vq, vi = 
(A 2 — 3A + l)vo. So if A is one of the roots in (|3.38|) . then there are solutions v spanned by v = (1, 1 — 
A, 0, some infinite word). Continuation of the iteration in (|3.36j) yields 

(1, 1-A, 0, A-l, -1, -1, A-l, 0, 1-A, 1, 1, 1-A, ...) 

(V Vi) (v 2 V 3 Vi V 5 , V 6 , Vr, V S , Vg),{ Vio, Ull),... 

Setting £ = (1, 1 — A) and rj — (0, A — 1, —1, —1, A — 1, 0, 1 — A, 1) the desired conclusion follows. □ 

These infinite-energy eigenvalues may be of physical significance as they serve to show that there is a 
natural way to "renormalize" , introducing a weighted sequence space, in such a way that in the renormalized 
Hilbert space, these infinite-energy eigenfunctions turn into finite-energy. 



4. Resistance metric 

Reviewing Definition 12.11 formula (|2.4p for the graph Laplacian, Ag, c the reader will note the coefficient 
c. Its significance is seen for example from electrical network models built on a given graph G = (G^°\ G^). 
Here c represents conductance, and it takes the form of a positive function defined on the set G^ of all 
edges. It signifies the reciprocal of resistance. The graph G is then a "large" (means infinite!) system of 
resistors, each edge e = (xy) representing a resistance of c(e) _1 Ohm between the neighboring vertices. The 
question arises of determining some sort of resistance metric giving the resistance between an arbitrary pair 
of vertices, so an arbitrary pair of points x and y in G^°\ This is non-trivial as there typically are many 
paths of edges connecting x with y. Nonetheless, the rules for electrical networks, Ohm's law combined with 
Kirchhoff 's law, allow us to compute a useful measure resistance, which we shall refer to as the resistance 
metric; see [Jor081 [Kig03l IPow76j . 

One of the conclusions in the operator approach to resistance amounts to identifying the resistance metric 
as a norm difference: The norm is now referring to a Hilbert space, and derived from the graph Laplacian 
Ag )C , i.e., the operator (|2.4|) for a particular choice of conductance function c. Hence there is a Hilbert space, 
the energy Hilbert space E = E(Ac, c ) and a function v from G^ into E such the resistance between x and 
y is the norm difference \\v(x) — v(y)\\E, i.e., the i?-norm difference between the vectors v(x) and v(y) in E. 
But each vector v(x) is itself a function on the set of vertices G^- ', in fact an electrical potential defined by 
Ag, c ] see Lemma l4~2l below. 

We further note that there is a general theory of metrics which can be computed this way with the 
use of Hilbert space, and we refer to the paper |Fug05| for an overview. In fact, these metric embeddings 
form a subclass of a wider family: spirals, or screw functions, as defined by von Neumann; and they have 
applications in information theory. In the present context, we need them for making precise the resistance 
metric, and to motivate the Energy Hilbert space. 
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Definition 4.1. Let G = (G^°\ G^) be a graph with conductance function c : G^ 1 ' — > R+. Pick a point o 
in G (0) . For x £ G (0) , the equation 

(4.1) &G,cV x = S - S x 

has a unique solution v x in the Hilbert space E obtained by completion relative to the quadratic form 

(4.2) H|:= 52c(xy)\u(x)-u(y)\ 2 . 

The expression £{u' ,u) is determined by polarization from (14. 2p , i.e., 

£(u',u)= E Yc{xy)(u'{x) -u'(y))(u(x) -u(y)). 

The resistance metric on G^ is given by 

dist(x,y) = \\v x -v y \\ E , (x,y € G (0) ). 
Here we are assuming that G is connected. For details see | Jor08] . 
Lemma 4.2. The solution v x in (|4.1| is determined uniquely up to an additive constant by the formula 



Proof. We check that 



if (v x ,u) = u(o) - u(x), (uEV). 



£(v, u)= 2 (A GyC v , u), 2 ( G (o)) 



Hence for all u G 2?, we have 



-^(i^m) = (Au x , u) /2 = ( by (@TT]) ) ((5 - S x , u) n = u(o) - u(x). 



□ 



Below we compute the resistance metric for the ./V-adic graph, and for the operator A GiC from Proposition 
EH and Theorem EHU 

Remark 4.3. If c : G^ — > R+ is a conductance function and if i> : G^ -* — > R is a function on the vertices 
of G, then (Ohm's law!) 

(4.3) I(xy) := c(xy)(u(x) - v(y)) 

defines a current flow on G. 

For a fixed pair of points x and y in G*- -' , we are interested in the following experiment which inserts one 
Amp at x and extracts it at y: It induces a current flow / : G^ — > R which is governed by (Kirchoff's laws) 

(4-4) J2 njk)=S x (j)-S v (j) 

Lemma 4.4. Let x,y € G^ ^ as above. If I : GW — > R is a solution to (|4.4p . and a voltage potential is 
determined from (|4.3[) . i/ien 

and 

(4.6) A G , c i> = (5 X - <*„. 

Proo/. First, (|4"3|) follows from a substitution of (|4.3p into (|4.2p . For l|4.6p . we get 

(a GiC tO(j) = E c 0' fc ) ~ ^( fc )) = ( by 63) ) E = ( by 63) ) = - s y (j). 

□ 
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Proposition 4.5. Consider the N-ary tree G = (G^ ', G^ 1 ') in Definition \3.1\ Let r\ := r\\ . , .rj n be a word 
in G?(°), n > 1. Then the solution (potential) v G E to 

(4-7) At> = <5 - <5^...^ 

is given by 

(4.8) l>(cJl . . . LU m ) = n — p{uJ\ . . . LU m , T]i . . . Tjn) 

where p(u>i . . . cu mi rji ■ ■ ■ T] m ) is the length of the largest common prefix for uo\ . . . io m and r\\ . . . r\ n , i.e., the 
largest p > such that p < m,n and u>i — r\{ for all i = 1, . . . ,p. 
The resistance metric is 

(4.9) dist(x,y) = y/2l(x,y), 
where l(x,y) is the length of the shortest path from x to y. 
Proof. We check (|4.7[) . 

N 

(Au)(0) = JVu(0) - = Nn ~ £ W W " viVi) = Nn - (N - l)n - (n - 1) = 1 = <J (0) - J,(0). 

i=l i^rn 

Next, we check the prefixes r\\ . . .r\ v with p < n — 1 

(Aw)(7?i ...ri P ) = (N + l)v(rfr...ri P )- . . .r} p i) - v(rfr . . .r] p+1 ) - v(rfr . . .7? p -i) = 

i¥=Vp+i 

(N + l)(n — p) — (N — l)(n — p) -(ra-p-1)- (n - p + 1) = = 5 $ (r)x . . . rj p ) - S„(jn . . 

(JV + 1) • - N ■ - 1 = -1 = 6 $ (rn . . . r) n ) ~ S v ( m . . . Vn ). 

All the other words have the form u> := r)\ . . . rj p uj p+ i . . . u m for some < p < n, oj p +i =f= Vp+i (if P = n ? 
then we just take m > n). Since the function v is constant n ~ p on the subtree with root r\\ . . . rjpWp+i, all 
the terms in the definition of (Av)(r/i . . . i] p ujp + i . . . u) m ) are equal n — p, and they sum to (N + 1) — N — 1 = 
= - 5 v (w). 

Clearly, only a finite number of edges (xy) have v{x) — v(y) ^ 0, namely the ones of the form rji . . . T) p ~ 
T)x . . . rjp+%. Therefore £(v, v) < oo. This proves (|4.8|) . 

To prove (|4.9p we use (Jpr08, Proposition 5.15]. Let p be the length of the longest common prefix of x 
and y, and let n, m be the lengths of x and y respectively. 

By |Jor08l Proposition 5.15] 

dist(x,y) = V2(v x (y) + v y (x) - v x {x) - v y (y))% , 

where v x is the solution (potential) for Av x = 5$ — 5 X and v y is the solution (potential) for Av y = 6$ — 6 y . 
Then 

dist(x, y) = V2(n - p + m- p- 0-0) 5 = V%l(x, y) 5 . 

□ 

The following corollaries involve two general issues for the resistance metric on graphs. They are motivated 
by two fundamental Hilbert space constructions: The first (von Neumann and Schoenberg) is a necessary 
and sufficient condition on a metric space (V, d) for the metric d to be the restriction of a Hilbert-norm 
difference. This then yields an isometric embedding of V into some Hilbert space. In the present case, this 
Hilbert space will be concrete: The energy Hilbert space E in Definition 14.11 and Lemma 14.21 The idea 
goes back to von Neumann and Schoenberg: The a priori condition on the metric d is that d 2 is negative 
semidefinite. The second construction amounts to Kolmogorov's consistency rules [PS72] , With this, we get 
a Gaussian stochastic process indexed by G^ . Kolmogorov [PS72] : An a priori given covariance function is 
positive semidefinite if and only if it comes from a stochastic Gaussian process. Below, we apply this to the 
energy-inner product (v x , v y ) E where for each x in 

G (o) 

v x is the potential function from Lemma [4~2l and 
(• , -) E is the inner product in the Hilbert space (|4.2p in Definition 14. II 
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Corollary 4.6. Let N e N and let G = (G^.G' 1 ') be the tree of Definitions^ For x,y e G<-°\ let l(x,y) 
denote the length of the shortest path from x to y. Then for all finitely supported functions £ S G^ - 1 — > C 
satisfying J2 X £{ x ) = 0, we have 

(4.io) EE^(*^(^° 

(In short the function l(, -) 2 is negative semi- definite.) 

Proof. Combining Lemma 14.41 and Proposition ^. 51 wc see that 

idist(a;,y) 2 = ~| 

where v : G (0) -> £ =(the energy Hilbert space), i>(x) := w x is the function determined in Lemma 14.21 
But it is known that the property in (|4. 10|) . (negative semi-definite) characterizes those metrics which have 
isometric embeddings into Hilbert space; see |Fug05| . □ 

Remark 4.7. The function v : -> E =(thc energy Hilbert space) from (|4.1ip also determines a system 
of covariances with the use of the inner products 

(V X , Vy) E = 2 (AV X , Vy} [2 = 2 (6 - S X , Vy) p = 2(Vy(o) ~ Vy (x) ) . 

An application of Kolmogorov's theorem |PS72j yields a Gaussian stochastic process w : G^ — > L 2 (Q,P) 
such that 



(4.11) l(x,y) = - dist(a:,y) 2 = -\\v x - v v f E 



(4.12) (w x ,w y ) E = / w x {w)w y {u))dP{w), (i,s6G (0) ); 

Jn 

and we may normalize with 

(4.13) / w x (u>) dP(w) = 0, (x G G (0) ). 



Corollary 4.8. Let N e N and let G = (G (0) , G (1) ) be the N-ary tree of DefinitionWIS. Then the stochastic 
process in (|4.12p has independent increments: specifically, if three points x, y, and z in G(°) satisfy x < y < z 
relative to the natural (partial) order of G^ ' then 

(4.14) £{v x -Vy,v y -v x ) =0. 

Proof. The order relation for the three points means that there are words r, s such that y = xr and z = ys, 
where we use concatenation of finite words. 

Using Lemma l4~2| we now see that (|4.14|) follows from the following identity 

v y (y) + v z (x) ~ v v (x) - v z {y) = 

which in turn follows from (|4.8p in Proposition 14.51 □ 



4.1. A random walk. In Theorem 14. 1 1 1 we will compute the moments of the measure /i c +p from Theorem 
13.261 For this we will use the following Proposition: 

Proposition 4.9. Let G = (G^°',G^) be a graph and c : G^ — > K + be a conductance function. We 
assume that G is not oriented, but it may have loops at some vertices, i.e., edges of the form (xx). Assume 
in addition that 

(4.15) ^c{xy) = l, (.xeG (0) ). 
Define the operator M. on T> by 

(4.16) (Mu)(x) = J2c(xy)u(y), (x€G (0 \ueV) 

yr^x 

Then A G , C = Ip - M. 

Define the random walk on G by assigning the transition from x to y, y ~ x the probability c(xy). For 
x,ye G<°> and n € No, let p(x, y; n) be the probability of transition from x to y in n steps. Then 



SPECTRAL THEORY FOR DISCRETE LAPACIANS 



23 



(4.17) M n s x = £ p(».y;»)*». (xeG {0) ). 

Let o € be a fixed vertex. Let fiM be the spectral measure associated to the operator M. and the vector 
5 Q . Then 

(4.18) J x n dfiM — (S , A4 n 8 ) — p(o. o; n) the probability of return to o after n steps. 
Proof. From the definition of M. we get that 

(4.19) M5 X = < x v) 6 ^ ( x e G(0) )- 

Equation (|4. 1T[) follows by induction: for n = 1, it is just the definition of M.5 X - Assume the equation is 
true for n, then 

M n+1 s x = Yp( x 'y'' n )^2 c (y z ) s z = Y^2 p ( x > y '' n ") c (y z ^ Sz = + 1 ) s ^- 

y z~y z y~z z 

Equation (|4.18p follows directly from (|4.17|) . □ 

Definition 4.10. Let T be the iV-ary tree from Definition 13.11 We define the graph T by considering the 
same vertices and adding an edge from to itself. 



Figure 2. The graph T (N = 2). 



00 ... 

01 ... 

10 ... 

11 ... 




Theorem 4.11. For each b£N, let Nj,(n) = the number of paths of length n inT from to itself. Then 

(4.20) J x n dfi c+p = —^=—Nf(n), (n e N). 

Proof. Let Ay be the Laplacian associated to the tree T with conductance constant 1 on every edge. Let 
AAf be the operator defined in equation (|4.16p for the graph T with conductance on each edge. Then 
a simple computation shows that 

(4.21) At = (N + 1)1 - (N + l)M f . 




From Theorem 13.261 we know that the At restricted to the cyclic subspace generated by Sq is unitarily 
equivalent to an operator of multiplication by N + 1 — 2^/Nx on L 2 (/i c+p ). Therefore the restriction of 
Aif to this cyclic subspace is unitarily equivalent to an operator of multiplication by ■ Using now 

Proposition ^. 91 we obtain that 



c+p — probability of return to after n steps, in the graph T. 



But this probability is equal to n^ro; and the result follows. □ 

Concluding remarks. We give the complete spectral picture for graph Laplacians for a number of 
classes of infinite graphs. Our conclusions include spectral representation, spectral measures, multiplicity 
tables, occurrence of semicircle laws, and rank-one perturbations. As corollaries of our main results, we 
note in particular that the spectrum must necessarily be absolutely continuous for the graph Laplacians Aq 
(with normalized conductance) when G is an infinite tree graph, and when G is a lattice graph. It would 
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be interesting to delimitate those infinite graphs G and associated conductance functions c for which the 
graph Laplacians Aq^ c have absolutely continuous spectrum. But it seems unlikely that a complete spectral 
representation is available for the most general infinite graph. 

A cl ass of graphs between two extremes is those generated by tree-state automatons. They are studied in 



GZ02) . and some spectral data is known. The class includes the free group Cayley graphs as well. 

We have explored two such graphs: (a) the graph of free group on two generators, and (b) those with 
vertices on Bethe lattices |Hug95| . Here the situation is more subtle: Case (a) is more complicated than 
the tree graphs in section 3 above. In section 3 we could build the spectral picture of A up with the use of 
a semicircle transform and a rank-one perturbation. But to get the spectral picture for A associated with 
the free groups, we might have to perturb by a rank-two operator, and the vector used for the perturbation 
might not be the obvious vacuum vector. 

The co-authors thank many colleagues for suggested improvements on an earlier version of this paper. 
The second named author (PJ) learned of these graph problems from Robert T. Powers in the 1970ties. The 
applications from that time include electrical networks of resistors, and KMS states in statistical mechanics. 
Since this, the interest in discrete analysis and graph Laplacians has grown tremendously. As a result, there is 
now an enormous literature, which in turn branches off in a variety of different directions. We cannot begin to 
do justice to this vast and diverse literature. It includes discrete Schrodinger operators in physics, information 
theory, potential theory, uses of the graphs in scaling-analysis of fractals (constructed from infinite graphs), 
probability and heat equations on infinite graphs, graph C*-algebras, groupoids, Perron-Frobenius-transfcr 
operators (used in models for the internet); multiscale theory, renormalization, and operator theory of 
boundaries of infinite graphs (current joint work between PJ and Erin Pearse.) We learned of some more 
such directions after the completion of this paper from groups of experts, and we are thankful to many 
colleagues who took the time to explain them to us. 
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man, Qiyu Sun, and Myung-Sin Song. The second named author thanks University of Central Florida for 
hospitality during a research visit where most of this work was done. Professors Strichartz and Golenia 
kindly helped us with literature citations after we completed the first version of our paper. 
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